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A first-principles derived approach is combined with the inverse Monte Carlo technique to deter- 
mine the atomic orderings leading to prefixed properties in Pb(Sco.5Nbo.s)03 perovskite alloy. We 
find that some arrangements between Sc and Nb atoms result in drastic changes with respect to the 
disordered material, including ground states of new symmetries, large enhancement of electrome- 
chanical responses, and considerable shift of the Curie temperature. We discuss the microscopic 
mechanisms responsible for these unusual effects. 



Perovskite alloys with the general form (A'A"...)(B'B"...)03 are of great technological interest because of their large 
piezoelectric and dielectric responses Examples of their applications are piezoelectric transducers and actuators, 
non-volatile ferroelectric memories, and microelectronic devices. Interestingly, little is known and understood about 
the microscopic mechanisms responsible for their very convenient properties. 

One particular aspect that seems very promising both to gain fundamental knowledge and to produce even more 
efficient materials is the effect of atomic ordering on the properties of these alloys [^-^). For instance, one may wonder 
whether some specific atomic arrangements can generate drastic and useful changes with respect to the usually-grown 
disordered material. Examples of such changes are a large enhancement of the electromechanical responses of the 
system or moving its Curie temperature into a temperature range of interest. If such arrangements exist, identifying the 
microscopic origin of their unusual properties would definitely be of great fundamental and technological importance. 
Addressing these issues in an efficient manner requires (1) a computational technique able to accurately determine 
the properties of sufficiently large alloy configurations and providing at the same time a deep microscopic insight; and 
(2) the ability to find the right atomic ordering leading to a given prefixed property, i.e. to solve the so-called inverse 
problem ||. 

The purpose of this Letter is to demonstrate that such an ideal situation is nowadays possible. More precisely, we 
have used the recently-proposed first-principles derived effective-Hamiltonian method of Ref. to tackle the inverse 
problem of designing and understanding ferroelectric alloys with prefixed properties. We focused on Pb(Sco. 5^0.5)03 
(PSN) alloy because Sc and Nb are chemically very different. (Note that Sc and Nb belong to different columns of 
the Periodic Table, which makes PSN an heterovalent alloy). As a result, we expected that atomic ordering in PSN 
might lead to significant changes with respect to the disordered case. As a matter of fact, we found that some 
arrangements between the Sc and Nb atoms (i) greatly enhance the electromechanical responses, (ii) lead to currently 
unobserved ground states of orthorhombic and monoclinic symmetries (while the disordered material adopts a well- 
known rhombohcdral ground state), and (iii) can considerably shift the Curie temperature. Furthermore, the analysis 
of our findings suggests that the mechanism responsible for these anomalous features in PSN can be understood by 
simple electrostatic considerations and/or simple phenomenological models. Some promising directions for future 
experimental and theoretical research are also provided. 

Modeling PSN- We first describe the effective-Hamiltonian approach of PSN. One starts by using the so-called 
"virtual crystal approximation" Q to create a simple system Pb(B)03 in which (B) is a virtual atom involving a kind 
of average between Sc and Nb atoms. The relevant lattice distortions of Pb(B)03 are described via a first-principles 
derived effective-Hamiltonian as the one used in Ref. |)| to simulate the finite-temperature properties of barium 
titanate. This Hamiltonian is then corrected to account for the presence of the two different B atoms in the real alloy. 
The complete Hamiltonian has the form: 

H eS ({ui}, {v,}, {77/}, {erj) = 

ffofftWJv^^D+ABHuJ.lv,},^}), 1 > 

where is the local soft mode in cell i, {v^} are vectors related to the inhomogeneous strains, and {rji} are the 
homogeneous strains in Voigt notation. The soft mode u is defined as the supercell average of the local modes: 
u = jfJ2i u ii where N is the number of cells in the system |T0] . The variables {cj} characterize the atomic 
configuration: <7j = +1 (resp. —1) indicates there is a Nb (resp. Sc) atom in cell j. 
Essentially, the correction term for PSN has the form: 
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where i runs over all the cells and j over the first three nearest-neighbor shells of cell i. Here, Bj% is the unit vector 
joining the B-site j to the B-site i; {Qji} are the alloy-correction coefficients, which only depend on the distance 
between sites % and j. Interestingly, Eq. (0) can be rewritten as the interactions between the dipole moments Z*Ui 
and radial fields Sji — the latter being generated by the tjj, being dependent on {Qji}, and acting on the neighboring 
cells i — : 



where £j can be viewed as the total field acting on u^. We find that the first-principles calculated Qji parameters all 
have negative signs. As a result, the field Sji created on site i by a Sc (resp. Nb) atom sitting on a site j is attractive 
(resp. repulsive). This is consistent with an electrostatic picture of PSN resulting from the fact that Sc +3 (resp. 
Nb +5 ) atoms are negatively (resp. positively) charged with respect to the average B-atom valence of +4; i.e. to some 
extent Sji can be viewed as the electric field created in site i by <jj. 

Equations (0)-(^) indicate that there is a competition in PSN between the global tendency of the system towards 
the ferroelectric ground state (GS) given by which is of rhombohedral symmetry with the local modes all 

oriented along the same (111) direction and the local tendency of the modes to align parallel to which typically 
favors an inhomogeneous local-mode distribution and therefore a decrease of the magnitude of u. Recent Monte Carlo 
simulations using this effective-Hamiltonian description of PSN have demonstrated its high accuracy to reproduce 
direct first-principles results in disordered, rocksalt-ordered (RS), and other non-trivially ordered structures 
These simulations also reproduce very well the finite-temperature experimental results for disordered PSN, once 
the simulation temperatures are linearly rescaled so that the calculated transition temperature coincides with the 
measured one PJTl|]. On the other hand, accurately predicting the properties of samples with local compositions 
significantly different from the 50% average of Sc and Nb may require the use of higher-order correction terms in 
Eq. (||). Even so, we believe that all the trends discussed in this Letter are qualitatively correct since Eq. (j^) includes 
the first-order (and thus the most important) correction terms. 

Solving inverse problems - The so-called inverse Monte Carlo is the most convenient computational scheme to 
search for the atomic ordering — i.e. the solution sample — leading to a prefixed property jq|. It is a Monte Carlo 
annealing done on the compositional variables Cj , and controlled by an energy functional 0[{ery}] that is defined so 
that its minimum — i.e. the result of the annealing procedure — corresponds to the solution sample we pursue. 
For instance, one would choose O = — | £^34 1 to determine the atomic configuration exhibiting the highest value of the 
(^34 piezoelectric coefficient. Note that we perform direct Monte Carlo simulations using the effective Hamiltonian of 
Eq. (Q) to evaluate the energy functional 0[{(7j}] in every inverse Monte Carlo step. Such evaluations usually require 
a long computing time and, as a result, we typically restricted ourselves to the use of a 4 x 4 x 4 supercell (320 atoms 
or 64 perovskite-like cells), which is smaller than the 12 x 12 x 12 supercell of Refs. HQ but still much larger than 
those feasible by direct first-principles methods. 

Tuning the GS structural properties of PSN - We estimate the GS structural properties (denoted in the following by 
the "0" superscript) by computing the corresponding thermal-averaged values at 5 K. The inverse problems considered 
here are listed in Table Q. We find that the RS sample maximizes the magnitude of the GS soft mode, which is consistent 
with the fact that the internal field £i is null in any cell i in this sample and, therefore, no inhomogeneous local-mode 
distribution resulting from Eq. (|j|) can exist. For all the other inverse problems indicated in Table [i], the solution 
sample can be described by a single non-null coefficient defined as: 



where k is a vector in the first Brillouin zone and Rj is the lattice vector associated to cell j |12[ . Our results show 
that the composition modulations described by k a = 1/4 = 3/4 (which correspond to the "Nb-Nb-Sc-Sc" sequence 
along the a cartesian direction) are related to the decrease of the a component vP a of the GS soft mode. This result is 
a direct consequence of the electrostatic-related Eq. (||) , since the other one-dimensional modulations compatible with 
the considered 4x4x4 supercell — i.e. k a = ( "Nb-Nb-Nb-Nb" or "Sc-Sc-Sc-Sc" ) and k a = 1/2 ("Nb-Sc-Nb-Sc") 
— result in zero internal field along the a direction and therefore in no reduction of u® a . Consequently, the sample 
defined by a single <7k with k = (1/4,1/4,1/4) exhibits the rhombohedral GS with the smallest magnitude of the 
soft mode. Similarly, the modulations described by k = (0,0,1/4) and k = (1/2,1/2,1/4) present a =><^" 
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sequence of very strong internal fields Si along the z direction. As a result, the z components of the local modes 
align parallel to these fields and the mean value u° z is null. This indicates that PSN samples with an orthorhombic 
GS can be produced by conveniently arranging the Sc and Nb atoms. We also tried to produce a tetragonal GS but 
only obtained a partial success, namely a sample characterized by k = (0, 1/4, 1/4) with a monoclinic GS for which 
3> \u^\ = \u®\. Interestingly, both orthorhombic and monoclinic phases have never been observed (yet) in PSN 
materials since the ground states of the two experimentally-grown samples (disordered and RS) are of rhombohedral 
symmetry Q. 

Maximum electromechanical responses at low temperature - We now consider the inverse problem of maximizing 



the electromechanical responses of PSN at 5 K 13 1. We focus on the piezoelectric coefficient and the dielectric 



susceptibility X33 because we numerically find that they exhibit much larger values than any other coefficient. In- 
terestingly, we also find that both of them reach their maximum values for the same atomic ordering. The solution 
sample exhibits %33 = 83 1 9 and ^34 = 4868 pC/N, to be compared with the values of 325 and 90 pC/N obtained in 
disordered PSN. A striking enhancement of the electromechanical response can thus be achieved by rearranging the 
Sc and Nb atoms. (We also obtained other samples with X33 > 4000 and ^34 > 2000 pC/N, all sharing the structural 
features that we discuss in the following.) The GS of the solution sample is characterized by u° = (0.544, 0.567, 0.076), 
where the main feature is the small value of the z component of the soft mode. This reflects the fact that the Ok 
coefficients that contribute most to the atomic ordering (around 70% of the total) have k z = 1/4. Our results thus 

suggest that samples with a "+ H " composition modulation in the z direction, where "+" (resp. "— ") stands 

for a plane that is rich in Nb (resp. Sc), will simultaneously exhibit small values of vP z and large values of X33 and 
c?34 at low temperatures. However, this modulation cannot be too strong since, for instance, we numerically found 
that the atomic configuration characterized by the single k = (0,0, 1/4) (for which vP z = 0, see Table Q) has small 
electromechanical responses. 

We now discuss the relationship between the small value of |u°| and the large electromechanical responses. Fig. |l| 
shows the probability distribution of u x and u z corresponding to the Monte Carlo simulation at 5 K of the solution 
sample, as well as the free-energy map of the system around u°. This latter is calculated out of the probability 
distribution as done in Rcf. fl4j| . The two main features of Fig. [j] are (1) the broadness of the u z distribution, which 
is responsible for the large value of X33 (see Ref. and (2) the flatness of the corresponding free-energy well, 

which indicates the ease of rotating the polarization around u° z . Interestingly, both features are typically associated 
with a small value of \u° z \ within the usual Devonshire-Landau description of ferroelectric perovskites p"5[ . Moreover, 
we numerically find that u z and 774 are strongly coupled in our effective Hamiltonian for PSN, which implies that 
their thermal fluctuations are very correlated. Consequently, according to Ref. Jl3|], a large value of X33 typically 
corresponds to a large value of CZ34, as we indeed find. 

Tuning the finite-temperature properties of PSN.- Studying PSN at finite (relatively-high) temperatures requires 
long Monte Carlo simulations in order to obtain good statistics. Thus, due to computational limitations we decided not 
to directly solve finite-temperature (FT) inverse problems. However, some links between the GS and FT properties 
are expected on the basis of well-known models that share some essential features with our effective-Hamiltonian 
description of PSN. More precisely, the so-called </> 4 model jl(| suggests the following relationship: the larger the 
value of \u |, the deeper the corresponding free-energy well at K and, therefore, the higher the Curie temperature 
T c . Hence, we thought it was pertinent to check whether the Curie temperature of PSN could be tuned by changing 
the value of the GS soft mode. Consequently, we studied the FT behavior of four samples that all undergo a single 
cubic to rhombohedral transition but exhibit different values for their GS soft modes. Our results, shown in Fig. |2j, 
confirm that it is indeed possible to move the Curie temperature to a temperature range of interest by designing alloys 
with a prefixed value of It is remarkable that relatively-small changes in the GS soft mode result in such large 
variations in T c . Interestingly, the transition of the sample with the smallest \u°\ takes place near room temperature. 
It is important to notice that one can tune the FT electromechanical responses of the system by moving T c , because 
these responses reach their maximum values around the paraelectric to ferroelectric transition. 

In summary, we have implemented a method to solve inverse problems in perovskite alloys by using the first- 
principles effective-Hamiltonian approach of Ref. M, and applied it to PSN. We predict that the system properties 
(GS symmetry, electromechanical responses, Curie temperature, etc.) can be greatly modified by varying the atomic 
ordering. Furthermore, all the anomalous features we predict can be essentially understood by simple electrostatic 
considerations and/or simple phcnomenological models. As a result, (1) we have indicated which kinds of alloy 
structures may lead to currently-unobserved structural phases as well as to giant dielectric and piezoelectric responses, 
and (2) we expect similar results to occur in other heterovalent perovskite alloys. Finally, we hope that our predictions 
will stimulate attempts at experimental growth of such samples. 
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TABLE I. Inverse problems related to several GS properties of PSN. First column: energy functional O. The "R" subscript 
indicates a rhombohedral symmetry (i.e. ur, = mr(1, 1, 1)). Second column: solution configurations. "RS" stands for the 
rocksalt structure while we indicate the k point leading to the non-null coefficient (see text) for the other cases. Third 
column: GS soft mode u°. In one case, we indicate two "good" solutions. We include the information about the disordered 
sample for comparison. 






sol. conf. 


u° (Bohr) 


-\<\ 2 


RS 


(0.507,0.507,0.507) 


\<\ 2 


(1/4, 1/4, 1/4) 


(0.431,0.431,0.431) 




(0,0,1/4) 


(0.552,0.552,0.000) 




(1/2,1/2,1/4) 


(0.484, 0.484, 0.000) 


K) 2 + {ulf 


(0, 1/4, 1/4) 


(0.665,0.242,0.242) 




disordered 


(0.461,0.461,0.461) 
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FIG. 1. Probability distribution of u x and u z at 5 K (solid line, in arbitrary units) and the corresponding free-energy map 
(dashed line, in ksT units indicated in the vertical axis) for the sample with the maximum responses discussed in the text. To 
obtain these results, we considered a 12 x 12 x 12 supercell (trivially obtained from the 4x4x4 solution sample by simple 
translation) and performed a 10 6 -sweep Monte Carlo simulation. Only the region around u° is explored, and the results for u y 
are essentially identical to those for u x . 
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FIG. 2. Temperature evolution of (mr) (see caption of Table |) for four PSN samples. These samples all undergo a single 
(cubic to rhombohedral) phase transition, and are chosen so as to exhibit different values of the GS soft mode. The values of 
it R considered here are: 0.507 Bohr (RS, which exhibits the largest possible value of UpJ, 0.484 Bohr (a intermediate sample 
generated by solving the appropriate inverse problem), 0.461 Bohr (disordered sample), and 0.431 Bohr (the sample defined by 
k = (1/4, l/4,a;l/4), which exhibits the lowest possible value of u%_). All the results correspond to the 2 x 10 5 -sweep Monte 
Carlo simulation of a 12 x 12 x 12 supercell, which thus guarantees a good convergence in the calculated transition temperatures. 
The temperatures have been rescaled so that the calculated T c of the disordered sample (944 K) coincides with the experimental 
value (373 K). 
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